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Abstract 

The quantum Maxwell theory at finite temperature at equilibrium is studied on compact 
and closed manifolds in both the functional integral- and Hamiltonian formalism. The aim 
is to shed some light onto the interrelation between the topology of the spatial background 
and the thermodynamic properties of the system. The quantization is not unique and gives 
rise to inequivalent quantum theories which are classified by #-vacua. Based on explicit 
parametrizations of the gauge orbit space in the functional integral approach and of the 
physical phase space in the canonical quantization scheme, the Gribov problem is resolved 
and the equivalence of both quantization schemes is elucidated. Using zeta-function 
regularization the free energy is determined and the effect of the topology of the spatial 
manifold on the vacuum energy and on the thermal gauge field excitations is clarified. 
The general results are then applied to a quantum Maxwell gas on a rt-dimensional torus 
providing explicit formulae for the main thermodynamic functions in the low- and high 
temperature regimes, respectively. 
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1 Introduction 

Quantum field theories at finite temperature have been intensively studied during the last years 
[U [2] . Prominent examples which stimulate this interest are the study of matter formation in 
the early stage of the universe, the description of the quark gluon plasma in the context of 
the AdS/CFT correspondence [3] and the current intense and partly controversial discussion 
concerning the Casimir effect [3] at finite temperature (for a review see El [HJ [9j ED]). 
From a general perspective, the thermal Casimir effect can be regarded as a deviation of the 
vacuum energy and the energy of thermal excitations of a quantum field caused by the presence 
of external constraints. These constraints may be imposed either by real material boundaries 
or by topologically non-trivial manifolds on which the quantum fields reside. In any case, the 
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modes of the quantum fields are correspondingly restricted, affecting both the vacuum energy 
and the thermodynamic functions of the system. 

In particular, it is the relation between the topology of the spatial background and the 
thermodynamic properties of gauge fields which in our opinion deserves closer attention and is 
the main motivation for the present paper. Our purpose is to study quantum Maxwell theory at 
finite temperature at equilibrium on a n-dimensional compact, closed and connected manifold 
X, which represents the spatial background. We perform the analysis in both the functional 
integral and Hamiltonian approach and derive the expression for the regularized free energy. 

Let us now motivate our intention of the present paper in more detail: Quantum fields at 
finite temperature on general manifolds with and without a boundary have been considered for 
many years PIIJ [13 Q31 021 [IEI [13 EDI Ell E21 E31 E31 E51 - (Here we have selected a list 

of references which are the most relevant ones in the context we are interested in) . Most of the 
research is in fact dedicated to the study of scalar fields. It is often argued, however, that the 
thermodynamic functions for a photon gas could be obtained directly from those for a massless 
scalar gas just by multiplication with the number of independent polarization states. We will 
show explicitly that this statement is not true in general. This requires a critical review of the 
concept which underlies the description of finite temperature gauge fields. 

In the functional integral formulation the thermal partition function of a bosonic systems 
at finite temperature T is obtained by integrating the classical action over all fields periodic in 
the (imaginary) time coordinate with period /3 := ■j^f, where ks is the Boltzmann constant 
[TJ[2]. Geometrically this corresponds to a quantum field theory of bosons on the product space 
To x X, where T 1 ^ denotes the circle (1-torus) with circumference /?. 

Evidently, the partition function plays a major role in the finite temperature context, since 
it serves as basic quantity from which the thermodynamic functions are derived. Hence par- 
ticular care has to be taken in order to obtain the correct measure for the functional integral 
representation of the partition function. This issue has been explicitly discussed in the seminal 
papers |26|, [27] for quantum Maxwell theory in the covariant gauge in the non-compact but 
topologically trivial case X = M. n . Following the Faddeev-Popov approach, the corresponding 
Faddeev-Popov determinant which appears as a field independent but temperature dependent 
multiplicative factor after gauge fixing must be retained in the partition function. Compared to 
the zero-temperature case, this factor is necessary to compensate the contributions caused by 
the unphysical degrees of freedom. 

But what happens if the thermal gauge theory suffers from Gribov ambiguities, which pre- 
vent the existence of a global and unique solution of the gauge-fixing condition, and thus the 
existence of any globally defined measure on the space of thermal gauge fields? 

It is a general result that abelian gauge theories suffer from the Gribov problem whenever 
the gauge fields reside on a non-simply connected, compact manifold |28[ 129]. However, this is 
precisely the situation one encounters in the finite temperature context since the base manifold 
is the product space T^xX. Thus even for topologically trivial spatial manifolds X, this problem 
does exist. This fact had been recognized long time ago |30[I31] in the case of quantum Maxwell 
theory at finite temperature on a 3-sphere (X = S 3 ) by solving the corresponding gauge-fixing 
condition. The main result was that the Faddeev-Popov operator (i.e. the zeroth order Laplace 
operator on x § 3 ) possesses temperature dependent zero-modes. Based on the results of 
Ref. |26j it was then argued, but not proved, that the Faddeev-Popov formula for the partition 
function is valid as long as the zero modes are ruled out from the domain of the Faddeev-Popov 
operator. Although this argument is valid for X = S 3 , we will show, however, that restriction to 
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the non-zero modes alone is not sufficient to give the correct partition function in the general 
case. 

In this paper the construction of the partition function and thus the Gribov problem will 
be tackled in a different way. As a by-product the relation between the occurrence of the zero 
modes and the Gribov ambiguities is elucidated, too. Thereby we will focus on the following 
four objectives: 

• Resolve the Gribov problem and construct a reasonable functional integral representation 
of the free energy including a finite expression for the vacuum energy using zeta function 
regularization. 

• Analyze the relation between the functional integral quantization and the Hamiltonian 
(canonical) quantization at finite temperature. 

• Study the impact of the topology and geometry of the spatial manifold X on the thermo- 
dynamic properties of the system. 

• Determine the thermodynamic functions for a photon gas confined to a n-dimensional 
torus (i.e. X = T n ). 

To our knowledge no comprehensive treatment of the various aspects of the Gribov problem in 
the finite temperature context has been given in detail so far. 

This paper is organized as follows: In Section 2 the geometrical structure of the space 
of thermal gauge fields is presented. It is shown that there exist gauge transformations not 
connected to unity. As a consequence the bundle of thermal gauge fields over the space of gauge 
inequivalent fields is not trivializable, which in physical terms is expressed by the statement 
that the theory suffers from Gribov ambiguities. 

Section 3 is devoted to the construction of a functional integral representation of the partition 
function in the space of thermal gauge fields. In order to circumvent the Gribov problem we 
will apply a method which has been introduced some time ago in the stochastic quantization 
scheme of Yang-Mills theory |32|, I33j. For abelian gauge theories this procedure was elaborated 
recently in Ref. [29], where it was also explicitly shown why the conventional Faddeev-Popov 
procedure fails on non-simp ly connected manifolds. 

The idea is to select a family of functional integral measures on the space of thermal gauge 
fields, whose domains of definition are determined by local gauge fixing submanifolds and which 
are integrable along the orbits of the gauge group. Finally these local measures are glued together 
in such a way that the physical relevant objects become independent of the chosen regulariza- 
tion of the gauge group and of the particular way this gluing was provided. The redundant 
gauge degrees of freedom are taken into account by factoring out the regularized volume of the 
full gauge group. This will be a slight extension of the original procedure [32], where only those 
gauge transformations were ruled out, which act freely on the space of thermal gauge fields. As 
a consequence of this reduction process, an additional topological factor apart from the con- 
ventional Faddeev-Popov determinant will appear in the functional integral measure. Whereas 
this factor may be neglected in the zero-temperature case, we will see that it contributes to 
the thermodynamical properties of the system and must be retained in the finite temperature 
context. Additionally, different topological sectors may exist on a general manifold X. In order 
to include all relevant thermal contributions, these sectors have to be considered as well when 
constructing the thermal partition function. 
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In accordance to the non-trivial topology of the space of inequivalent thermal gauge fields, 
the quantization is not unique. In the functional integral approach this fact is usually taken into 
account by adding a total derivative to the classical Maxwell Lagrangian. This term does not 
alter the classical equations of motion but contributes in the quantized version whenever the 
integration runs over topologically non-trivial field configurations. In our case this additional 
action will be shown to be parametrized by a constant vector £ R^W, where bi(X) is the first 
Betti number of X. The functional integral can be solved exactly and by using zeta function 
regularization technique [Ml E3 ES] we obtain a closed expression for the free energy of the 
system. This is achieved by expressing the zeta-function of the Laplace operators on x X in 
terms of the zeta function associated with the corresponding Laplace operators on X. 

The quantum Maxwell theory at finite temperature is discussed from the canonical (Hamil- 
tonian) point of view in Section 4. We will determine the true phase space V and provide an 
explicit parametrization of this manifold. Since V turns out to be non-simply connected, there 
exist inequivalent quantum theories which are classified by unitary irreducible representations 
of tti(V). These are labelled by the non- integer values 6 £ M bl ^ x \ We will determine the Hamil- 
ton operator and calculate its spectrum. The finite total vacuum energy of the system turns 
out to be the sum of the vacuum energy of the transverse gauge fields, which is regularized 
by adopting the minimal subtraction scheme [35], and the energy of the ^-states. The latter 
vanishes whenever 9 G Z bl ^ x \ The free energy is derived within the Hamiltonian scheme and 
the agreement with the functional integral formalism is explicitly proved. 

Finally, the high-temperature asymptotic expansion is calculated in terms of the heat kernel 
(Seeley) coefficients giving rise to finite size and topological contributions to the familiar Stefan- 
Boltzmann black-body radiation law. 

The scaling property of the system under constant scale transformations is analyzed in 
Section 5. Due to the regularization ambiguities of the infinite vacuum energy and the non- 
trivial cohomology of X, the free energy does no longer transform homogenously. This leads to 
a modified equation of state. 

In Section 6 we discuss the example of finite temperature quantum Maxwell theory on 
the n-torus X = T n . The low- and high temperature expansions are calculated for the main 
thermodynamic functions in terms of the Epstein zeta function [37, 38J and the Riemann Theta 
function. Due to the #-term the ground state is degenerate. In the zero-temperature limit the 
entropy converges to the logarithm of the degree of degeneracy, which proves explicitly the 
validity of the 3rd law of thermodynamics (i.e. Nernst theorem). 

Necessary results on Riemann Theta functions and Epstein zeta functions are summarized 
in the Appendix. The conventions c = h = ks = 1 are used. 

2 The configuration space of thermal gauge fields 

In this section we will consider the geometrical structure of the space of thermal gauge fields. 
These results are based on Ref. [29], where a detailed analysis can be found. 

Geometrically, the finite temperature quantum Maxwell theory at equilibrium at finite tem- 
perature 1//3 is regarded as gauge theory on the product manifold Tjxl with product metric 
9 = 7/3 © 7- Here denotes the 1-torus T 1 equipped with the temperature dependent Rieman- 
nian metric 7^ = /3 2 dt <g> dt (t is the local coordinate on T 1 ) and the spatial background X is a 
n-dimensional compact, connected, oriented and closed Riemannian manifold equipped with a 
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fixed metric 7. The corresponding volume forms on and X are denoted by vol T i and volx, 
respectively. 

Let (Q, ttq, x X) be a principal £7 (l)-bundle over T^xl with projection 7Tq. The space of 
thermal gauge fields is identified with the C°°-Hilbert manifold A® of all connections of Q. The 
space of M-valued fc-forms on x X, denoted by Q k (Tj 3 x X), is equipped with the L 2 -inner 
product 

<vi,v 2 > g = v l A* g v 2 v 1 ,v 2 e n k (Tl x X), (2.1) 

JtlxX 

where * g is the Hodge star operator on x X, satisfying i? g = (— l) fc ( n + 1_fc ) when acting on 
it-forms. The co-differential d\ = (-1)(™+ I )( fc + 1 )+ 1 * g d n+1 ^ k * g : ^(T^ x X) O fc_1 ( T J x x ) 

T 1 xX 

gives rise to the Laplacian operator A fc s = d* k+1 dk + d^-id^,. Let % fc (Ti x X) denote the 
space of M- valued harmonic A;- forms on x X and 1-L k (T^ x X) -1- its orthogonal complement. 
Then we can define the Green's operator by 



<V lO^xll^^xI) 1 , :=(A/ |^ (T i x^)" 1 ^ (T ^ xX) , (2.2) 

where n v 1 3 ' is the projection onto 77. (T^ xX) . By construction one immediately obtains 

A f\xX TlxX TlxX f),xX W lt/ T l „ yU 1, 1 

A fc 8 o G/ = G fc 8 o A/ = n H ( V A) . The dimension of % fe (T^ x X) is given by 
the k-th Betti number 6fc(T^ x X) = 6fc(X) + bk-i(X). Here bk(X) is the fe-th Betti number of 
X. 

Let t 1 := V — 1 K denote the Lie-algebra of £7(1). The space of thermal gauge fields .4.^ is 
an affine space modelled on fi x (T^ xljgt 1 and can be equipped with the following metric 

g(wi,W2) =< -j=, -j= > g (2.3) 

where ^1,^2 £ ^(^ x X) t . Let vol^ Q denote the induced (formal) volume form on 

The gauge group Q is defined as the group of vertical bundle automorphisms on Q and 
can be identified with the Hilbert Lie-Group C°°(T^ x X; £7(1)) of differentiable maps. The 
corresponding Lie algebra is <5 = Or(Tg x X;! 1 ). Since there is at least one non-contractible 
loop in x X, the gauge group contains gauge transformations which cannot be connected to 
the unity. These transformations are classified by a non-trivial ttq(Q) and cannot be generated 
by taking the exponential of imaginary- valued functions on x X. Under an arbitrary gauge 
transformation v E Q, the thermal gauge fields transform according to 

A ^ A v = A + (tt* q v)*# u W v£ g, (2.4) 

where £ ri 1 (T 1 ;t 1 ) is the Maurer Cartan form on £7(1). (For notational convenience 

we shall not distinguish between tt*qV and v.) The isotropy group of the non-free action f|2.4[) 
is £7(1), which corresponds to the subgroup of constant gauge transformations. However, the 
quotient group := Q/U(l) provides a free action onto AP giving rise to a smooth manifold 
= AP /G* which represents the reduced configuration space in the functional integral 
picture. 
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In order to fix notation, let %\{X) denote the abelian group of harmonic R- valued differential 
/c-forms with integer periods along Zfc(X;Z), where Zfc(X;Z) is the subcomplex of all closed 
smooth singular /c-cycles on X with integer coefficients. 

In our case we can choose a set of 1-cycles a G Z\(X, Z), i = 1, . ..,6i(X), whose corre- 
sponding homology classes [cj] provide a Betti basis of ffi(X;Z). Let p< n-1) G ^^(X) be 
a basis associated to the Betti basis [q] via the Poincare duality, where j = 1, . . . , 6 n _i(X). 
A dual basis G %\{X) can be adjusted such that J* C .P£ = fx Pi ^ P^f ^ = Let 
(h!x)ij =< Pi,Pj >7= fx Pi ^ *7 /4 denote the induced metric on "^(X). Notice that 
*iPi = Y^j 1 =i \^x)ijPj l • The harmonic I-forms (pr^ pr x p^\ . . . , pr x p^ x ^) gener- 
ate U\(J\ x X), where = ^ wI T i. Here pr T i : T 1 x X -> T 1 and pr x : T 1 x X -»■ X are 
the natural projections. 

Let (sq,xo) G x X be an arbitrary but fixed point. The gauge group G* is equivalently 
characterized as particular subgroup {v G G\v(sq,xq) = 1} of Q. The corresponding Lie algebra 
<S* = (3/t 1 contains all those £ G (5 such that £(so,xo) = 0. We use this alternative description 
to rewrite any (restricted) gauge transformation as product of an infinitesimal- and large gauge 
transformation. For this we construct an isomorphism k : (5* x H^Tp x X) — > by 

a)(s, x) = exp £(s, x) • exp (27iV— 1 / «)• (2-5) 

Here C( s : [0, 1] — > x X is an arbitrary path in x X connecting (so,xo) with (s,x). As 
a result, (|2.4[) can be rewritten into the form 



A i-> A + d £ + 2vr\/^T I m pr*i + ^ m jP r x py> j , £ G ©*, m , G Z. (2.6) 

For any choice of an arbitrary but fixed background gauge field Aq G A® there exists a smooth 
surjective map vr^°g : M? -> T 1+bl W, defined by 

*£ ([4) = ( e /T ^^> (A - Ao \ e / {so}xci (^o)___ ) / {so >x CblW (^o)^ (2?) 

According to the general result proved in [29], we thus obtain the following two propositions, 
which summarize the topological structure of the space of thermal gauge fields: 

Proposition 1. (A®, tt^q , M$) is a non trivializable flat principal Q^-bundle over M.® with 
projection tt^q . □ 

Proposition 2. For an arbitrary but fixed connection Aq G AP , '■ -M.* — > T 1+&1 ^ admits 
the structure of a trivializable vector bundle over T 1+bl ( x ) with projection tt^Jq o- n d typical fiber 

M := imd\ (git 1 . If A' Q G .A^ is a different background connection, then : Af? —5- T 1+6l ( x ) 
is an isomorphic vector bundle. □ 

As a result, the manifolds A® and T 1+bl W x N x are locally diffeomorphic. An explicit 
expression for these local diffeomorphisms can be given as follows: Let us introduce the two 
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contractible open sets V aj= i = T l \{northernpole} and V aj= 2 = f l \{southernpole}, which 
cover the j-th 1-torus T 1 C T 1+bl W. Then 



V = {V a := V ao x ••• x V aj x ••• x V abi(x) \ a := (a , . . . , a j: . . . , a bl{x) ), ^€{1,2}}, (2.8) 

is an open cover of Hence the family of open sets U A ° := { , k j ^q)^ 1 {V a ) provides a 

finite open cover IA A ° = {U A °} of M*. Any two local sections a a . : V aj — > R 1 of the universal 
covering bundle R 1 T 1 generate 2 1+bl W local sections a a : V a C J 1+b ^ x ) K 1+il ^ , given 
by a a = (cao , ■ ■ ■ , Ca bl {x) ) • It can be shown that the maps Xa° '■ Va x AA x — > tt^q (U a ° ) C A® , 
defined by 

/ bi{X) \ 

Xa°(^o,---,^ l( x),T,v) = A) + 2ttV z 1 U ao (zd)pr;ie+ ^ ^(^i^pfM + r + u*^W. 

(2.9) 

provide an appropriate family of local diffeomorphisms. 

3 The free energy in the functional integral approach 
3.1 The construction principle 

In this section we want to construct a functional integral representation of the (thermal) par- 
tition function for the photon gas. The Maxwell fields are governed by the classical action 

S mv (A) = \\\F A \\ 2 = 1 -g(F A ,F A ). (3.1) 

In the previous section it has been shown that the reduced configuration space is the non-simply 
connected manifold A4?, since 7Ti(A4?) = kq{G*) = It is well known that ambiguities 

arise in the definition of the vacuum if the corresponding configuration space is non-simply 
connected. In the functional integral approach this fact is usually taken into account by adding 
a so-called topological- or theta action, denoted by Sg, to the gauge invariant classical action. 

To each arbitrary but fixed vector 9 = (0i, . . . , &bi(x)) £ K^ 1 ^-' we associate the harmonic 
1-form := ^n^^^jSi^x )!/ 6 '^ e ^(X) t 1 . We propose the following action 

Se(A) ■ = 9(FA,* g pr* x * 7 0) 

bi(x) (3.2) 

i=i J l\* x 

which is indeed topological since it is independent of the metric of x X. This specific choice 
will become more transparent in the next section when the Hamiltonian approach is considered. 
For the total action we take now 5 tot := Si nv + Sg instead of Si nv alone. Since ^ = 0, the 
inclusion of this additional term does not change the classical equations of motion. The theta 
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action Sg is imaginary because our context is the finite-temperature Euclidean regime where 
time is compactified. In fact, analytic continuation of (|3.2p onto the real time axis is necessary to 
obtain the correct real- valued theta-term in the Hamiltonian formulation on space-time Rxl. 
The total configuration space, denoted by Aji xX , is the disjoint union 

A nxX = \jA® (3.3) 
Q 

over equivalence classes of principal [/(l)-bundles Q over the base manifold x X. Each class 
is labelled by its first Chern-class c x (Q) G H 2 ^ x X;Z) H^X;^) H 2 (X;Z). 

In the conventional functional integral formulation the integration of the functional 
e~ Stot vol 9 A Q over A® would become infinite due to the gauge invariance of the classical ac- 
tion. Our aim is to give this integral a well-defined meaning by damping the contributions from 

the gauge dependent degrees of freedom. The non-trivial bundle structure A® aQ ) MS allows 
only for a local damping procedure, which, however, can be patched together using a partition 
of unity in the end. 

T 1 xX 

The map [£] \-> d*doG ^ (£) provides an isomorphism (5* = imd* <X> t . Notice that Q is a 
trivializable principal C/(l)-bundle over whose global trivialization is given by the diffeomor- 
phism a(u) = (uu(so, xo)~ l ,u(so,xo)) where u G Q. Let -d G G ^(G) <8> <3 and G ft 1 ^*) <5* 
denote the Maurer Cartan forms on Q and G* , respectively. Formally one can introduce the nat- 
ural metrics T g = g($ G (.), $ G {.)) and T g * = g($ e * (.), $ G * (.)) on these two gauge groups. Let 
volg and volg t denote the induced left-invariant volume forms. Furthermore U(l) is equipped 
with the standard flat metric giving rise to the volume form uoZj/m = {\/ "& u ^\ With 
respect to our parametrization of U(l), the corresponding volume is Jj/m voljj^ = 2n. 

Let us now introduce three real- valued regularizing functions A G eg G C°°(Q), A G * g G C°°{G*) 

and ArjP G C°°(U(1)) in such a way that the following correspondingly regularized group 
volumes 



Vol(Q;e A ^ e ») := / volg e < oo 
Jg 



Vol(Q*;e A " e * 9 ) := / volg^ e < oo (3.4) 

Jg* 

A u(1 \ f , -A u{1) 



Vol(U{l);e- Are « ):= / vol u(1) e' 1 ^ 
Ju(i) 

become finite. Evidently, a regularization of U(l) is not necessary (of course Areg = would 
be sufficient) but we will use this freedom later on to absorb irrelevant multiplicative factors in 
the partition function. 

In order to relate these volumes, we consider the differential of <t -1 , which gives 



Td--'(Z,w) = Tr^(Z + -\ t= ove w J , feTA, w G T Z U(1), (3.5) 

where r G denotes the right multiplication on Q by z G U(l). Since = i?^*(£) G irad\ (g) t 1 

the induced metrics are related by 
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(£-i)*T° = T e * - pv $ u{l \.)$ u{1) {.), (3.6) 

where V := J x volx is the volume of X with respect to the metric 7. For the volume forms, 
this finally implies 

(a^Yvolg = (/3V) 5 uolg, A uoZtf (1) . (3.7) 

Let us take the natural choice A^ eg := TTgKf* g +a1 x ^A^eg \ where irg : Q — > G* is the projection 
and o"( SO]a;o )(w) := u(so,xo)- Using f|3.4j) it follows that 

Vol(g-e- A ™9) = (pV)* Vo/(^;e- A ^)Vo/([/(l);e- A ^ 1) ). (3.8) 
An appropriate choice for the regularizing function is provided by 

Afe» = ^||^K^ (1) )H 2 + ^\\U nl ^ xX \v*^)\\ 2 , (3.9) 

where LT^ ( T /9 X ^) is the projector onto 7^ 1 (T^ x X) [29] . According to f|2.5|) and with respect 
to the chosen Betti-basis, each v 6 £/* is uniquely represented by a pair (£, m) G (imd* <8> t 1 ) x 
ji+bi(X) ^ ^ s a conse q Uence the integration over G* splits into an integral over imd\ (g) t 1 and a 
summation over the components of the vector m = (mo, • • • , fnbi{X)) G Z 1+6l (^. For the choice 
(|3.9p one obtains 

Vo/(^;e- A ^) = (detA^^l^o^x^x)-^^^)^^^^^^), (3.10) 

where Sf )1 (x)(-|-) denotes the 6i(X)-dimensional Riemann-Theta function (]A.lj) . 

Now we come to the construction of the functional integral representation for the thermal 
partition function. Due to the non-trivial bundle structure of the space of thermal gauge fields 
the separation in gauge independent and gauge dependent degrees of freedom can be done only 
locally. The idea is to construct integrable measures locally in A and in the end to paste them 
together with a partition of unity to obtain a global and integrable measures. This method 
was originally introduced for studying the relation between stochastic quantization and the 
conventional Faddeev-Popov quantization scheme [32] and then applied to define an integrable 
partition function for Yang-Mills theory in order to overcome the Gribov problem [33]. Recently 
this approach was used for the formulation of an appropriate functional integral representation 
of generalized p-form gauge fields [39J . 

Let {p a } be a partition of unity of M.* subordinate to U A ° and define oj^° := prg t o ip^°, 

where ip£° : 7r^g(C/ a ) — > U a X G* is a family of local trivializations of the bundle A® aQ > M.® 
and prg t : x G* G* is the natural projection. Let us introduce the following smooth 
functional on A®, 

a 

which accounts for the regularized group volume G*- Now we define the functional integral 
representation for the thermal partition function by 
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Z((3,V) = N V f V ° 1aQ s ~ Q e~ St °K (3.12) 

The partition function is the formal sum over all equivalence classes of principal U(l) bundles 
over Tjxl classified by their first Chern-classes c\{Q) and the functional integration over the 
corresponding spaces A® of thermal gauge fields. The normalization constant N is temperature- 
and volume independent and will be fixed later on. The volume of the total gauge group Q is 
factored out in order to reduce the system to the true physical degrees of freedom. Notice that 
this is a slight generalization of the original procedure introduced in Ref. [33], where only the 
volume of the restricted gauge group was factored out. It will be shown that this proposed 
partition function is integrable and resolves the Gribov problem. Moreover, the phase e~ Se gives 
the weight-factor for the different topological sectors. 

Physical observables are regarded as gauge invariant functions on Aji xX - The thermal 

expectation value (TEV) of a physical observable / G C°°(A®) is defined by 

Ec(Q)gH2( T i xA:;Z) Ifl(/) 

< f >p= wx) ' (3 - 13) 

where 

if (/) = f "f*\ - e~ s ^ E Q f . (3.14) 



p Jaq Voi(g- e - A ^) 

As a consequence the TEV of a physical observable is independent of the particular choices of 
the regularization, the local trivializations and the partition of unity (for an explicit proof see 



3.2 Determination of the free energy 

In accordance with the bundle structure of A®, the functional integral in (|3.12p can be trans- 
formed into an integral over 1-1+6! (X) x N x g„ which can be explicitly calculated. Using the 
family of local diffeomorphisms (|2.9p . one obtains for the transformed measure 



(Xa°T vol A Q = ( deth ji x x det A o' 3XX |-H (TixX)-L ) uoZ T i+6 1 (x) \v a A volj^f A volg^ , (3.15) 



where volji+t^x) \v a is the induced volume form on f 1+bl ^ x \ yet restricted to the patch V a . 
It follows that volji+biix) = i^volu^ A ... A iXi(X) vo ^U{\) where ik- U(l) ^+ T 1+bl ^) denotes 
the k-th inclusion for k = 0, . . . , b\(X). The restriction of g to N induces a flat metric on 
that space and an associated volume form volj^. Apart from the conventional Faddeev-Popov 
determinant of the scalar Laplacian, the factor det h2i v appears in addition in the induced 

measure (j3. 15|) . This factor carries the topological information of x X and will be shown to 
depend on temperature and volume. On the contrary to the zero-temperature case, where this 
factor as well as the conventional Faddeev-Popov determinant may be neglected, it is essential 
to keep both factors in the functional integral in the finite temperature context. 
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A direct calculation gives 



1 



*s [prjiQ^'j = jj pr* x vol x , 



(3.16) 



so that the induced metric on T-L l {T\ x X) admits the following matrix of rank 1 + b\(X), 



,(i) P~ l v 
V* ~ { o phg 



Hence 



det h, 



(i) 

JlxX 



^iW-Vdet^. 



(3.17) 



(3.18) 



Due to the Hodge decomposition theorem it is always possible to choose the fixed background 
gauge field Aq G A® in such a way that the corresponding field strength Fa becomes harmonic, 
i.e. ^—j=Fa G %i(T^ x A). Since the cohomology of X is finitely generated, the Chern-class 
c\(Q) admits the following (non-canonical) decomposition 



'bi(X) b 2 (X) _^ 

ci(Q)= I £ £ m j » 7 f ) + £ rf } I eH 1 (X;Z)(BH 2 (X;Z), 



(3.19) 



i=i 



j'=i 



k=l 



with respect to the Betti bases ^ and {rfpfff for il^AjZ) and H 2 (X;Z), respec- 
tively. Here 62(A) = dim H 2 (X; R). The torsion subgroup TorH 2 (X; Z) is generated by the basis 
(ij. )fc = i- Hence there exists n, . . . , G N such that r^u = for each k = 1, . . . , w . The order 
of the torsion subgroup is \TorH 2 (X; Z)| = rjLi r k- Finally the vectors (Zi, . . . , l bl (x)) G Z 6l W, 
(mi, . . . , m^m) G Z b2 W and ?/£ G Z rfc = Z/r&Z denote the components with respect to these 
different generators. 

Let (p$ 2) )$ 2 jf } be generators for then (pr^ Apr^pf^/ and {pr^pf)^ are 

the induced generators for "Hf (T^ x A). The background field strength can then be expressed 
in the form 



F Ao = 2ttV^I 
A direct calculation gives 



'&1PO b 2 {x) 

i=l j=\ 



(2) 
.7 



(3.20) 



*g (pr*x pf) = P pr^Aprx (* 7 pf ) 

* s (p4^Wa^p«))=1 



(3.21) 
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yielding the following metric (h^ xX ) on / H 2 (T 1 ^ x X), namely 

V* - ^ o /34 2 V' (3 ' 22) 

of rank &i(X) + bi{X). Here (h^)ij =< p[ 2 \ >-y is the induced metric on T-L 2 (X). 

Let us recall that on an arbitrary compact manifold M the Laplace operator, when restricted 
to ^(Af)- 1 , splits into the sum A^f\ nk ( M )x = ^f\imd k ^ + ^f\imd* k+1 - Since the spectra of 
A^|j mc i fc _ 1 and A^Jimd* coincide, the spectrum of Aff\ H k(M\x is the union of eigenvalues of 
A^f\i m d* +1 and of A^ 1 |j mc ;*. In our case M = x X. This implies for the determinant 

fi x jf T 1 xX T 1 xX 

detA/ | w i (T i xX) ± = (det A/ \ imd *) (det A Q | w o (T i xX) x). (3.23) 

Using ([321), (13TT5L (EHHJ), (I3T20]1 . (I3T22D and (ET23]) . the integration in (13T41 can be carried out 
and gives (for / = 1) 



i9d) 



=(2tt) 



W) (detail 



(det A,/ 



■H°(TixX)- L )(det A x l^i(TixX) J 



x exp 



(27T) 5 



6i(X) 6 2 (X) 



+ 27TA 



_«*(*) 

-1 2 0^ 



x ^oZ(Z7(l); e - A - e » 3 )- x iV. 



(3.24) 



The pair of vectors (I, in) G x Z &2 W labels the inequivalent principal c7(l)-bundles Q. 

Using the decomposition (|3.19p . the sum over the Chern-classes c\(Q) is performed by summing 
(|3.24p over the free part and the torsion part of -ff 2 (T^ x X;Z), respectively. 



Let us choose e 



(2tt 



for the regularizing functional on U(l), then the thermal 



partition function ()3. 12[) admits the following form 



Z e (P,V) =(2vr)^/3^ (detfc«)i (det A? XX | w o (T i xX) x) (det A? X \ 1(T i xX) x)^ 
x 8 6l(x) (e|27TA/=T /T 1 /#>) e ta(x) (0|2vrV^T /3 hf) \TorH 2 (X;Z)\N. 

(3.25) 

In order to give the formal determinants arising in (13. 25ft a mathematical meaning, the zeta- 
function method will be used as regularization technique. Generally, the regularized determinant 
of a non-negative, self-adjoint elliptic operator T> of second order on a general manifold M is 
defined by 

detP = exp(-C / (0;P)) =exp^-^| s=0 C(s;P)). (3.26) 
Here C(s;^) denotes the zeta-function of the operator V, defined by 
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*.feo r(s) Jo 

where the sum runs over the non vanishing eigenvalues u a (T>) of T> only. In this sum each 
eigenvalue appears the same number of times as its multiplicity. This sum is convergent for 
9i(s) > dl ™ — . The second equation in (|3.27|) is the heat-kernel representation of the zeta- 
function based on the Mellin transformation. Therein UP is the projector onto the kernel of T>. 
Let V denote the restriction of V to the non-zero modes, i.e. V := 'D\ ket . v ±, then £(s;P) = 
C(s;T>') holds by construction. In general, the operator D possesses the following asymptotic 
expansion 

oo 

+T> v ^ — "S . , k — dim M 

Tr(e~ w ) ~ J2a k (V) t — a — , for U 0, (3.28) 
fc=0 

with respect to the asymptotic sequence t i— > t 2 of functions [40\ ETT] 142] . The constants 
a k (T>) are the Seeley coefficients of T>. By using this expansion and splitting the integral in 
(|3.27p into an integral over [0, 1] and [1, 00) respectively, the zeta-function admits the asymptotic 
expansion 



/ 00 \ 

Eafc(2?) , a diTaM {V) -dimkerP 

- 1 fc-dimM + " + r(s;D) 

\k^dimM ) 



(3.29) 



where r(s;D) is an analytic function. The zeta function has a meromorphic extension over C 
with simple poles at s k = dim for k G No and residue Res s=Sk [£(s; V)} = at 

s = Sfc. Since lim^o ^(s) = 1, the zeta function is analytic in the origin and one finds that 
C(0;P) = adimM(^) — dimkerP. In our case this leads to 

C(0; A? XX ) = C(0; A? X \ P(T i xX)± ) = a n+1 (A? XX ) - b p (T} x X) for p = 0, 1. (3.30) 

Let us denote the eigenvalues of A^ by i/^(A^), where l s runs over an appropriate subset 
j( s ) C Z. Notice that the eigenvalues of A r ' 3 (r = 0,1) are given by u k {A r p ) = (^p) 2 with 

T 1 xX 

k £ Z. Then the spectrum of Ap 13 \w(T 1 x.X)- L (P = 0> 1) contains the following set of non- 
vanishing real numbers 



bpec{A p \ n p(Tl:_ 

(3.31) 

:= v k {K P ) + ^(Af ) ^ 0|r + s = p, (fc, Z s ) G Z x ■/«}. 



In terms of these eigenvalues, the corresponding zeta functions can be written as convergent 

2 



infinite series (for 9J(s) > ^s— ) 
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C(*;A?* X 



Z)(fc,Z )GZx J<°) 

Z](fc,« )eZx j(°) 



,(0,0) 
(Mo) 
,(1,0) 

(Mo) 



(Jc,ii)eZx j(!) 



,(0,1) 

(Ml) 



(0,0) 



,(1,0) 
(Mo) (k,loY 



if p = 0, 

ifp = l, 

(3.32) 



where each eigenvalue appears as often as its multiplicity. Remark that v { 

We notice that there is an intrinsic ambiguity in the definition of the zeta function (|3.27p due 
to the fact that the eigenvalues of the Laplace operators are not dimensionless. As a consequence 
the partition function would admit a dimension as well. To restore this, an appropriate scale 
factor fj, of mass dimension [p] = 1 has to be introduced [33]. Instead of (|3.27p one introduces 
therefore the dimensionless zeta function £p(s;D) := J2u a (v)^o(t I ~ 2l ' a ('^ > ))~ S ■ Formally, this 
gives (^(s;V) = // 2s £(s;P). In the following all determinants appearing in ()3.25|) are expressed 
in terms of these scale dependent but dimensionless zeta functions, i.e. 



det (A/ 



HP{TjxXy 



exp 



— \s=o CtM{s; A p p 



p = 0,l. (3.33) 



However, even with this substitution the partition function f|3.25j) is not dimensionless. This 
traces back to the occurrence of the two Jacobian determinants (|3.7p and (|3.18p . whose dimen- 
sions must be corrected accordingly. In our convention the coupling constant is set to 1, which 
leads to the mass dimensions [A] = and [Fa] = for the components of A and Fa, 
respectively. Consequently, one obtains the mass dimensions [h^ ] = —1 and [h^] = 1. 

In order to get a dimensionless partition function we propose the following replacements in 
(1377]) and ([57lTD . namely 



pv ' y f3Vfj, n+1 

h$ xX 1 ^ diag (/T 1 VfjT- 1 ,^//) 



(3.34) 



Alternatively, this substitution could be implemented formally by choosing the normalization 
constant N = //iW- 1 in ([335]) . 

As a result, the free energy J*(/3;V) = - 4 In Z e {/3;V) admits now the form 



2/3 
1 

2^ 



1 

1 



InG 



61 (X) 



0|2tta 



r 1 fcS } ) 



-^lne 62(x) (0|2vrV^l /3 



In 2?r - 4 In \TorH 2 (X: Z) 

2/3 /3 1 v ' 7 



1 , . TixX, s 1 _ JixX, . l-&i(X)\ 

+ |-C(0;V | w o (T i xX) x)-— C(0;A/ Ui (T i xX ^) + ^^)ln/^ 2 - 



/3 



2/3 



(3.35) 

In the next step we want to express the free energy in terms of the geometry of the spatial 
manifold X: Given any non- negative, self-adjoint elliptic operator D of second order on X with 
spectrum {u a (D)\a 6 J}, we define the following two auxiliary quantities 
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(k,a)e(ZxJ)' 



l(s;D): 



keZaeJ' 



(3.36) 



The prime in the first formula indicates that the sum runs only over those indices (k, a) such 
that (^|p) 2 + v a (D) 7^ 0. In the second line the index set J' Q J labels the eigenvalues of the 
restricted operator D' . One immediately finds that 



l(s;D) =l(s;D) + 2(dimkerD) 



2tt 



2.s 



Cfl(2«), 



where Cr( s ) is the Riemann zeta function. Hence f|3.32[) can be rewritten in the form 



C(s;A 



p \w{r^xxy 



'l(s;A x ) ifp = 0, 

X( S ;A^)+X( S ;Af) ifp = l. 



(3.37) 



(3.38) 



Applying the Mellin transform to X(s; D) and performing the integration over t give 



X( S;j D) 



r(*) 



v a {D')t 



2^ T{s) 



(3.39) 



Here K u (s) is the modified Bessel function of the second kind [33]. In order to calculate the 
derivative of I'(0;D) in s = 0, we notice that C(s;.D) admits a Laurent series expansion in 
s = — namely 



C(S"2^) 



(3.40) 



fc=i 



where FP denotes the finite part of the zeta function in s = — ^. Recall that for an arbitrary 
meromorphic function /, the finite part in sp S C is defined by 



FP s =s [f] := lim 



/(-) 



ds. 



t->0+ 27TV-1 J\ s - S0 \=e S~ S 

(see e.g. |44| for the definition and properties). For the zeta function of D this implies 



FP s=S0 [((s;D)} =]im - 

e— >Q Z 



C{-\-e-D) + a- l - + e ] D) 



(3.41) 



(3.42) 
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Taking the expansions Jjv = s + 7s 2 + 0(s 3 ) and T(s - 5) = r(-|)(l + *(-|)s + 0(s 2 )) for 
small s, where 7 is the Euler-Mascheroni number and is the Digamma function |43j . and using 
that the relation lim^o Js(f^j) = h(0) holds for any regular function h, one finally obtains 



l'(Q;D) = - 2(dimker J D)ln/3 - 2 ^ In [l 



e 

aeJ' 



- B FP s= _i Ms; D)] + 2(1 — In 2)Res, [((s; D)] 



(3.43) 



T 1 

Using the duality formula ()A.3|) the heat kernel expansion of A p fi becomes 

rr(e --?) = e l( 0|V=I^) = ^=6,(01^) -^+0(.-*), forUO, (3.44) 
for p = 0, 1. This implies afc(A p ' 3 ) = f^^fe.o- By comparing the asymptotic expansions of 

T 1 xX , . 

A p p and using (|3.3ip . one finds the following relation between the corresponding Seeley 
coefficients, 



a k (A ? ) = a k (A$) 
a k {A 1 } XX ) = J^[a k (Af) + a k {AX) 



(3.45) 



We want to calculate the sum in the second term of (I3.43|) in the case of D = Af. Let 
{i/ ai (Af \imd*)\ ai G Ji} and {z/ ao (A(f | imd *)| a G J } be the spectra of Af | w * and Afi\ imd *, 
respectively. Since the spectrum of A^ \^ P ^± is the union of eigenvalues of A^ \i m d* +1 an d 
Ap^limd* the sum runs over v ai (Af\ imd * 2 ) and v ao (A^\ imd *), respectively. When inserting the 
explicit expression (13.38P into (13.35P only the sum over the eigenvalues of Af \i m d* survives. 
Furthermore, this also implies ((s; Af \i m d*) = C{s;Af) — ((s;Aq ). In terms of the Seeley 
coefficients, the residue of £(s; D) in s = — \ is given by 

Res s= _,_ {((s; D)} = _ a <"+^) , (3 .46) 

2 2-y/7T 

Due to the linear structure of FP, Res and applying the duality formula (|A.3P once again to the 
6i(X)-dimensional Riemann Theta function, we finally get the functional integral representation 
for the free energy 



T'ip-.v) =\fp„^ [C(»:Af \ imds )] + \Re>,=-l [C(»; Af \ imd .)] ln(^ 



-/3^ Q1 (Af | imd .) 







-lne b2(x) [0|2~s 



-1/J(fcg>) 



1 

-^ine 6l(x) 

^ln|TorF 2 (X ; : 



2tt 



(3.47) 
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It will be shown in the next section that the Hamiltonian formalism gives exactly the same 
result. Eq. (j3.47p shows explicitly that the ambiguity in the free energy is independent of the 
temperature. 

Before closing this section we want to annotate briefly why the conventional Faddeev-Popov 
procedure is not applicable in the present case (For a detailed account refer to [29J). In the 
Faddev-Popov approach the starting point would be the functional integral (|3.12p . however 
with the choices A^, eg = A~!* g = 0. If the gauge fields were constrained by the covariant gauge 
condition d*(A — Aq) = 0, this would lead to the well known additional gauge fixing term 
^||<f[(yl — ^4o)|| 2 i n the action. This term in combination with the classical Maxwell action 

T 1 xX 

would then give the unrestricted kinetic Laplace operator acting on the thermal gauge 

fields in the total action functional. Since 'H 1 (T^ x X) ^ this operator is not invertible, so that 
the functional integral over the gauge fields would diverge. This is precisely the Gribov problem 
in the abelian case. Even if this problem was solved ex post by restricting the Laplace operator 
to the non-harmonic forms, the Faddeev-Popov determinant related to the covariant gauge 

T 1 xX 

condition would be - following the conventional approach - only the factor det Aq 13 \fi°(T^xX)- 
In our approach there is an additional multiplicative factor det hji x which turns out to be 

essential in the finite temperature context and which appears naturally as part of the Jacobian 
of the transformation (|2.9p . It gives rise to the first term on the right hand side of f|3.35|) . Its 
first part, namely jjln/3, is related to T-L l (T^) and is present even if b\(X) = 0. However, this 
term is compensated by the zero-mode subtraction coming from the expansion ()3.37p of the 
zeta function of the Laplace operator on x X in terms of the zeta function associated to 
the corresponding Laplace operators on X. Without this cancellation, the entropy (|5.ip would 
have admitted a logarithmic divergence in the zero-temperature limit . In this context we would 
like to refer to the scalar case, where this logarithmic term was present in |17] but disappeared 
in |23j due to a modification of the functional integral formula for the free energy. By this 
modification the zero modes of the scalar Laplacian are taken into account correctly and the 
equality with the operator (Hamiltonian) approach is provided. 

The second contribution, namely fel 2 ^P ln/3 (i.e. the second term in (|3.35p ) is completely 
absorbed in the Riemann Theta function due to its duality property. Moreover, our treatment 
of the Gribov problem automatically rules out the zero modes of the kinetic Laplace operator 
right from the beginning, so that the functional integral can be carried out, yet giving a finite 
result. 

4 The free energy in the Hamiltonian approach 

In this section the quantization of Maxwell theory is studied from the canonical (Hamiltonian) 
point of view. The aim is to determine the corresponding Hamilton operator H 8 in the presence 
of #-vacua and to calculate the (thermal) partition function and the free energy according to 




(4.1) 



-lnZ e ((3;V). 
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The trace is taken along the physical (i.e. gauge invariant) states. If the theory possesses different 
topological sectors (like it is in our case), one has to perform in addition a sum over these sectors 
as well. We will return to this topic below. In a first step the objects in (|4.ip are marked with a 
caret in order to distinguish them from the partition function and free energy obtained in the 
functional integral formalism. 

As in the previous section, let Q be a principal [/(l)-bundle over x X and consider 
the pull-back bundle Q := i* x Q over A, where ix '■ X <-)■ x X is the canonical inclusion 
lx{x) '■= (1, cc). For the fixed-time canonical formalism we consider - as usual - the principal 
t7(l)-bundle R x Q over the space-time manifold R x X. 

Let A— denote the space of connections of Q and let A s = f2°(A; t 1 ) be the space of scalar 
gauge potentials. The tangent bundle T(A9- x A s ) = A^- x A s x ^(Ajt 1 ) x ^(A^t 1 ) is the 
corresponding configuration space, which is parametrized by the coordinates (A, A S ,A, A ). 
After performing the analytic continuation, the topological action Sq (I3.2D becomes real. The 
Lagrangian associated to the classical Maxwell action Stot admits the form 

L 9 (A,A S ,A,A S ) = l\\A - d A% - + ^7(i- d A s ,6), (4.2) 

where Fa = d\A is the magnetic field. Here |||| 7 refers to the norm which is induced by the 
metric 7(^1,^2) := - J* x t?i A* 7 t>2 on A^-, where vi,v 2 G 1 (A;t 1 ). Since ^(A-dAf^O) =A /(A,0) 
holds, the topological term is a total time derivative, so that the classical equations of motion 
remain unchanged despite the addition of this particular term. By performing the Legendre 
transformation one gets the corresponding Hamiltonian in the phase space T*(AQ x A s ) with 
conjugate momenta 

U= 6 S=A-dA s + -^, ir = ^J = 0. (4.3) 
- SA ~ ~ 2tt' - 5A S K ' 

The phase space T*(A— x A s ) is equipped with the canonical symplectic form. Let A £ A— 
be a fixed background connection satisfying d^FAg = 0. The non-vanishing Poisson brackets for 
the basic linear phase space functionals A u := j(A — A , u), A^ := j(A s , v), Tl u i := 7(n, v!) and 
II*/ := 7(n s ,T;') are given by 

{A U ,U U/ } = j( U ,u'), {Aim,} = j(v,v'), (4.4) 

where u,u' € fi 1 (A;t 1 ) and v,v' € ri°(A;t 1 ). The Lagrangian L e is singular and leads to two 
first class constraints in the phase space, namely IP 5 = and the Gauss law d*Tl = 0. By 
choosing the temporal gauge A s = 0, the first constraint is solved. The dynamical system is 
restricted to the submanifold C := {(A, II) G T*^^|<i*n = 0} and governed by the induced 
Hamiltonian 



H e (A,m = \\\U-^\\* + l\\FA\\ 2 T (4.5) 



The corresponding group of gauge symmetries is q0- = C°°(X; U(l)). Let Q— '■= G—/U(l) be 
the restricted gauge group, which acts freely on A—. In fact, the Gauss law constraint is related 
to the symmetry under the identity component G^~ of and can be solved by factoring out 
this subgroup. This leads to the quotient space C/Q^~ . However, since Q^r is not connected, 
i.e. ttq(G^-) = Z bl ( x \ this is not the true physical phase space V. That space is obtained by 
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taking out also the large gauge transformations. As a result V = (C /G~ )/tto(G~) = C/G~, 
which is symplectically equivalent to the cotangent bundle T*Mt of the space of gauge orbits 
Mir = A—jGiT- At first sight this seems to be not the final result, since only the restricted 
gauge group Q— instead of G— has been considered so far. But since 

A Q/qQ ^ M ^ T*Mt can 
be really viewed as the true physical phase space V of the classical system. 

In the Hamiltonian approach the underlying geometrical structure is the principal t/jr-bundle 
A— —7- Ai^r. In order to obtain an explicit parametrization of Aij-, we can apply propositions 
1 and 2 of section 2 (by replacing the base manifold x X with X in the assumptions). 

Thus Mt admits the structure of a trivializable vector bundle over T bl ^ with typical fiber 
Af := imd\ <8> t 1 , where d\: Q 2 (X;t 1 ) — > 1 (X;t 1 ). The trivialization induces a vector bundle 

isomorphism <p: T*(T bl ^ x ^ x Af) — > T*A4t between the corresponding cotangent bundles. If 
(zi, . . . , z bl ( X y,Pi, ■ ■ ■ ,Pbipf)) are coordinates for T*T bl ( x ) and (r, T) are those for T*Af, then 
4> is given by 



4>{zi, ■ ■ -,Zbi(x),T;pi, ■ ■ ■ ,p bl (x),Y) = 

bi(X) 



2vr 



T bi(X) 



j,k=i 



where the brackets denote the equivalence class in At— and s a . : V^. C T 1 — > R are the local 
sections of the universal covering R — > T 1 (as defined in Section 2). It can be easily shown that 
(j) is globally well defined. Notice that we have parametrized T*Ai~ by pairs (L4],n) satisfying 
d*n = 0. The inverse map reads 



cuu-jj.) (eix(^o)Apr i \... )e ix(^o)AP^; )d * G2(i ^_^ )) 

2 _ 6l r n_ M*) 



E("S , ) 1 */ x a^r , .-...^ T D'S'Im^/^a,'-", (4.7) 



<£JG 2 diIIJ . 

The field strength i 7 ^ = (L4 of the background gauge potential classifies the topologically 
non-trivial monopole configurations. With respect to the chosen Betti-basis (see section 2) one 
gets 



b 2 (X) 



F Ao = 2ttV^1 Y m kPk\ rn k G Z. 



(4.8) 



fc=i 



In these adapted coordinates the classical Hamiltonian splits into two independent dynamical 
subsystems with the phases spaces T*T bl ^ x "> and T*A[_, respectively. The dynamics is governed 
by the Hamiltonian 4>*H e = H e harm + H trans , where 
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H harm( z l' • • • > z bi(X)iPU • • • >Pbi(X)) = 

= E (w - **0i)iPi - ^j) + V- E )« ( 4 - 9 ) 

and 

#tran S (T;T) = 1 1| T || * + i 7 (r, Ai r). (4.10) 

In order to quantize these two subsystems we determine the Poisson brackets between the conju- 
gate variables. Let us define the angles qj = 27T lrr[ fx(— ~ — o)^P^ ^ f° r 3 = 1> • • • > &i(-^Qi which 
can be regarded as coordinates of (i.e. of the universal cover of T^W). By definition 

Zj = e 27Tyr ~^ qj and from (|4.4p we find {<2j,Pfc} = ^ifc- Each gauge transformation v G £/— of the 
gauge fields, i.e. A 1— >• .A*, induces a translation ^ \-¥ qj + ay, where ay = 27r ^- f I c V *$ U ^ £ ^ 
are the winding numbers related to v and j = 1, . . . , b\(X). The Poisson bracket of the linear 
phase space functionals t u := 7(7", u) and T u / := 7(Y,it') yields {r u ,T M /} = 7(u, g^G^iV)) 
where u,u' G il^Xjt 1 ). 

In the Schrodinger representation the Hilbert space of physical states splits into the 
tensor product S) = Sjharm &> fttrans of Hilbert spaces of sections of line bundles over T bl( - X ^ 
and over M_, respectively. The C/^invariant metric 7 induces a natural connection in the bundle 
A— — > Aij- by declaring the orthogonal complement to the fibers as horizontal subbundle. The 
metric restricted to this subbundle finally induces a metric 7 on A4~. With respect to the 

diffeomorphism M^r = T bl ^ x J\[_ this metric splits into the direct sum 7 = © lj\f, where 
1/V is the induced flat metric on M_. Let volj^ denote the induced volume form, then the (formal) 
inner product in S) is given by 



<*i,* 2 >i5= / voL bl ( X ) (dethy) 3 ipiip2 / voltf ipnp2, *i = A® <Pu i = 1,2. 

(4.11) 

The Hamilton operator fl 6 harrn is obtained from (14.9)) by substituting the classical momenta by 
the operators pi = - 7 =^-. This leads to 

-, M*) 2 b 2 (X) 
Harm = E ( fc SV + E ( 4 ' 12 ) 

where Vj := — 2tt\/— 19i is the covariant derivative. This system can be interpreted as 
quantum theory of a classical particle in b\{X) dimensions coupled to a constant (functional) 
electric field 9 and moving in an external constant potential determined by the topological 
sector m G Z b2( - X \ 

The energy eigenstates of H harm are the wave functions 4>j{qi, ■ ■ ■ , gwjn) = e ^ -^=1 i<?i 
with eigenvalues 
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Lrn 



+ 



(27T) 



2 62 W 



(4.13) 



fc</ 



In the previous section we introduced the topological action Sq and argued that this term 
accounts for the inequivalent quantum theories caused by the topologically non-trivial config- 
uration space. Now we will justify this argument: Let us introduce the unitary operator Uq 
by 



${q) := (U ij>M 



illi ■ ■ • ,Qbi(X))- 



(4.14) 



Since UoVjU a 1 = -Mr, the ^-dependent term in (I4.12p can be removed giving rise to a new 



den 



U e^harm U e 1 with eigenstates ipj{q) := Ugipf Since H e harm and 



Hamilton operator Hh arm 
Hharm have the same spectrum, the corresponding quantum theories are equivalent. However, 

the wave functions ip are no longer single- valued, since ip{q + a) = e ^j =1 aj j i()(q) 

for a £ 7j bl ( x \ For each fixed 6 E R bl W these wave functions can be regarded as sections 
of the line-bundle C e := x^-C over T bl ( x \ which is associated to the universal covering 

jj&iPO — y j; bl ( x } via the unitary irreducible representation a \-t e 27r v-ia0 Thus different choices 
for 9 ^ lead to inequivalent quantum theories. 

Now we are going to quantize the transversal modes: Let {ro^Jai G Ji} denote an or- 
thonormal basis of eigenforms satisfying A x \ im( i* zu ai = u ai (A x \i m d%) ro ail where each eigen- 
value ^(A^lj^*) appears as often as its multiplicity. With respect to the decompositions 
T = Sa ie Ji T «i CT «i and T = EcueJi T "i ro ai one obtains 



trans 



aiSJi 



IT 



tti 



r> fax (^1 \imdX 



' cti I 



(4.15) 



which is nothing but the (well-known) Hamiltonian of an infinite number of harmonic oscillators 
with frequencies J v ax (A x \i m d* )• The coefficients are r ai =< A— Ao,m ai >, T Ql =< T,m ai >, 



r ai and T ai 



having the properties r ai 

{T'aiJal} = 
commutators [r ai , T a / 
representation by introducing annihilation and creation operators 



T ai . The non- vanishing Poisson brackets are 



5 Q1)0 ,' i , which in the quantum theory are replaced by the 
1 fiai,a' ■ The spectrum of (I4.15P is easily obtained in the Fock 



J ai ■- 



(4^ ai (Af 



imd. 



s))' 



(4i/ ai (Af | imd * 



^qi (A^ limc/j) ^"ai 



(A^ | jmdj ) 



T 

- T 



(4.16) 



which have the non- vanishing commutators b ai , b a , = S ai )Q / . The transverse Hamilton oper- 
ator admits then the following familiar form 



i 

aiGJi 



(4.17) 
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which admits the energy eigenvalues s\^ s = yu ai (A^ : |j m <f*) (k ai + |), where k ai £ No- 

In order to determine the thermal partition function, we have to take the different topological 
sectors into account. These are labelled by the first Chern class (see (|3.19p ). namely 

b 2 (X) w 

ci(Q) = i\ Cl {Q) = £ mjrjf + ^ y k tf G H 2 (X; Z). (4.18) 

j=l k=l 

Since H^ arm depends on the free part of H 2 (X; Z) only, one obtains 



Z e (J3;V) = Tr { e 

C1 {Q)£H 2 {X;Z) 



E E II E e fce-^" \TorH\X;Z)\. 



(4.19) 



where the trace has been calculated by summing over all physical eigenstates. The free energy 
is then given by 



P{P;V) = \ £ ^MXrW + I E ln 1-e'V^^I^) 



--ln9 6l(x) 
- ^ln|Tor^ 2 (X;Z)|. 

Obviously, the first term in (|4.20p , which represents the vacuum energy € vac of the transverse 
modes of the electromagnetic field is infinite and requires a regularization. We choose zeta- 
function regularization and introduce the following function in s 6 C 

e(s) :=-T fc^l*"^)* = I^ + i c( Af j } (421) 
2^ (A - 2l/Qi(AfM)r+ | 2^ 

where /i is a scale parameter with mass dimension [/2] = 1. This step is necessary to assign e the 
correct mass dimension. A natural choice for the vacuum energy could be to take e vac = e(— \). 
However, this would be reasonable unless ((s; Af \i m d* ) becomes divergent in s = —\. In order 
to account for this case as well, we follow [35J and define the regularized vacuum energy of the 
transverse modes as finite part of e, i.e. 

z r :i--=FP s= _r\e(s\ (4.22) 



2 



This regularization scheme amounts to remove the pole from (|4.2ip . Since FP S=S0 [/i(s)/a(s)] = 
fi(so) F P s=SQ [f2(s)] + f{(so)Res s=so [f2(s)] holds for a function f\ being holomorphic in so and 
fi being meromorphic with simple pole at the same point sq, one finally gets 

<7c = \FP s= _,Us;AX\ imd *)] + l -Res s= _ h \as; A* U*)] ln£ 2 , (4.23) 
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which represents the finite vacuum free energy of the transverse modes of the Maxwell field. 
Let us now substitute the first term in (|4.20p by s V a C and choose fl = then one finds 

f e ((3;V) = F e {l3;V). (4.24) 

Hence we have explicitly verified that both quantization schemes are equivalent. Let us stress 
the fact that keeping the field independent Jacobian, taking the quotient by the volume of the 
total gauge group and finally summing over the different topological sectors have been the main 
steps in the functional integral scheme to obtain this equality with the Hamiltonian approach. 
In the remainder of this paper the caret on the free energy will thus be omitted. 

The normalization scale n expresses the ambiguity of the free energy which occurs whenever 
Res s= _i [C(s; Af-\ im( i* )] 7^ 0. In that case renormalization issues have to be considered. Alter- 
natively, this could be stated as follows: If £(s; A^|j m d*) is finite for s = — g, the free energy is 
uniquely determined and the scale dependency must disappear. 

In the case of massless scalar fields at finite temperature the relation between these two 
quantization schemes and in particular the question regarding the zero modes of the kinetic 
operator have been discussed some time ago in Refs. [181 US GO] and recently in Ref. [23J . 

In order to split T V) into a temperature independent term F^iV) and a temperature 
dependent contribution Tf e rnp (f3;V), we use (|A.lj) to rewrite 



(oi^(^r') - 

According to the modular property of the Riemann Theta function (|A.2p . 9 can be replaced by 
the translated vector (6) whose j-th component is defined by (6j) := sgn(6j) mm mj£ z \irij — 
where j = 1, . . . , b\(X). Hence the components of (0) are restricted to \ (6j)\ < | for all j. 

The expression for the free energy of the quantum Maxwell field at finite temperature admits 
now its final form 



6 



T t> (P;V) = T°(V)+T° emp (l3 ] V) 

= \ FP s=-l [Cfa A? U rf *)] + \Res s= _i [(is; Af \ imd .)] In (^) 
1 bl{X) 1 

j,k=l aiSJi 

- ilnG 62(x) (O^V^T^) " ^ln\TorH 2 (X;Z)\. 

This formula for the regularized (total) free energy of the photon gas confined to a closed 
manifold X is the main result of the present paper and - to the best of our knowledge - has not 



-/y^ Q1 (Af \ md *) 



(4.26) 
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been stated before. It exhibits clearly how the topology of X affects both the vacuum energy and 
the thermodynamic structure of the system. The vacuum energy J 7 ® is the sum of the regularized 
vacuum energy of the transverse modes and the ground-state energy of Hfl arm (see (|4.12p ). The 
latter vanishes whenever 6 G Z bl ( x K In any case the Laplace operators appearing in (I4.26h are 
correspondingly restricted in order to rule out any zero- modes. This is the consequence of the 
construction of the partition function by using either a family of local trivializations in the 
functional integral approach or an appropriate parametrization of the true phase space in the 
Hamiltonian scheme. 

The free energy is unique only if the zeta function converges at s = — A. On the other hand, 
if the first and second cohomology group of X vanish, the Riemann Theta functions as well 
as the temperature independent ^-vacuum term disappear and the free energy is completely 
determined by the transverse modes. Under these conditions the free energy of the quantum 
Maxwell theory would be indeed a multiple of the free energy of a massless scalar gas. The 
n-sphere X = S n with n ^ 1, 2 is a typical example for such a configuration. In so far we gave 
a proof for the statement argued in [30], [31] • 

In the case n = 1, the transverse modes are absent so that the free energy is exclusively 
governed by the harmonic component. For X = § 2 the ^-states are absent, but since H 2 (§> 2 ; Z) = 
Z, the topological sectors contribute additionally to the thermal excitations. 

Let us now compare our result with the free energy thermal contributions of a photon gas in 
flat Euclidean space confined to a very large box in Euclidean space. In fact, when considering 
the thermal excitations in the infinite volume limit, the index ot\ in the fourth term in the 
second equation in (14.26P becomes the continuous n-dimensional wave vector k. Hence all finite 
size and topological effects are neglected. Furthermore, the sum is replaced by an integration 
with respect to the measure . Using [33] one gets finally for the thermal part of the free 
energy density in that limit 



(n -1)J *L In (l _ e -*ffl) = _(„ _ T( 2+l) Cfi(n + l) ^-(n+D, (4.27) 

where the factor n — 1 is the number of independent degrees of freedom. Since the topology of 
Euclidean space is trivial the Riemann Theta functions are absent. Eq. (|4.27p is the well-known 
Stefan-Boltzmann term of black-body radiation. 

Let us emphasize that in our study the photon gas is confined to a closed spatial manifold 
X. Evidently, there is no space "outside" of X, which could contribute neither to the vacuum 
energy nor to the thermodynamic excitations of the system. In that sense, (I4.26P represents the 
intrinsic free energy of the photon gas. This has to be distinguished from configurations where a 
partition separates an inside and an outside region in an ambient space. The resulting (Casimir) 
force exerted by the photon gas on that partition is then caused by the difference - the so-called 
Casimir free energy - between the free energy in the inside and outside region. Typical examples 
are infinite parallel plates [10] , rectangular cavities [231 HI] an d piston geometries [35j UHJ H7] . 

As was stated above, in our case the vacuum (Casimir) energy depends, in general, on the 
normalization scale /x. However, if the electromagnetic field is confined to a compact and con- 
nected cavity with smooth perfectly conducting boundary within M 3 , the corresponding Casimir 
energy was shown to be finite [H]. This is based on an explicit computation of the heat kernel 
expansion of the corresponding Laplace operators and the fact that relevant contributions from 
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the inside and the outside of the cavity cancel. Hence for that configuration no renormalization 
is necessary. 

In the remainder of this section we want to determine the high-temperature limit of the free 
energy. The starting expression will be (|3.35p together with (|3.38[) . Let us define the auxiliary 
quantity 



K(s;D') :=2Y,Y, 

k=l a£j' 



2irk 



(4.28) 



where the second sum runs over all eigenvalues of D' , thus excluding the zero modes of D. By 
separating the zero-modes, the quantity I(s;D) in (|3.36p can be alternatively rewritten in the 
form 



l(s;D) = as;D') + 2(dimkerD) (J-\ ( R (2s) + K(s; D') 
The Mellin transformation of K(s; D') reads 

2 



K(s-D') 



T(s) 



Jo agJ 



(4.29) 



(4.30) 



k=l 



for Sft(s) > ^ but can be analytically continued elsewhere [19]. Its asymptotic expansion for 
f3 — > can be obtained by substituting the heat kernel expansion (|3.28|) for the restricted 
operator D' with corresponding coefficients, denoted by a m (D'). The integration term by term 
gives 



K(s;D)c2^ T j ^a m (D)^j) — . (4.31) 

In order to perform the limit s — > of K(s; D') and of its derivative respectively, we notice 
that the function T(s + n )(n(2s + m — n) has simple poles at m = n and m = n + 1. 
Using that E\(s; 1) = 2(jj(2s), where E\ is the Epstein zeta function (|B.1|) in one dimension, 
the reflection formula (IB.5P provides the analytic continuation of that function. If we take the 
Laurent expansion of (r(s) at the pole s = — 1 and the series expansion for (see above), a 
lengthy calculation yields 



--\ s=0 K(s;D')= a m(D') 2 n - m ir-h m - n T( n ~ m )C R (n + 1 - m) 

m=0 

"f™ (4.32) 



+ a n+1 ( J D')^r(ln(^)+ 7 ) -a n (D')ln/3, 



where (once again) 7 is the Euler constant. From a m (D) = a m (D') + 5 mjn (dimker D) and the 
spectrum of (p = 0, 1), it follows that 

a m (Af \ imd *) = a m (Af ) - a fc (A^) + (1 - b x {X))5 m , n . (4.33) 
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Using (fOll . (021) . (Q3D . the explicit form a (A^) = (4tt)~% , ( ^i p) , y [40J and finally apply- 
ing the duality formula (1A.3|) to the Riemann Theta function ©b 2 (jq(0|...), one obtains from 
(|3.35p the high-temperature asymptotic expansion for the free energy in terms of the coefficients 
(Af \imd*), namely 



J*(fi; V) ~ - (n - 1) 7T- 2 * 1 r(^±A) Oi(n + 1) /T (n+1) V 



- £ a m (Af | w *) 2™- m vr-2 / 3 



m— n— 1 



„ . n + 1 — m , , , 

r( r ) 0?(n + l-m) 



1 

2^ 
1 



C'(0;Af U*) + ln f ^^^ iTori^z/ 



(4.34) 



+ — [2a n (Af \ imd *) + h(X) + b 2 (X)] ln/3 



+ Res s= _i [C( S ;Af| iWS )] (ln(f^) + 



/3 



iine b2W (oi^r 1 ^)" 1 



The first term highlights the familiar Stefan Boltzmann term in n spatial dimensions (see 
(|4.27p ) and the remaining terms represent the modifications caused by the topology of X. If 
C(s; Af\i m d*) is regular at s = — \, the //-dependent term vanishes as expected and yields an 
unambiguous result. A corresponding formula for a gas of massless scalar fields on a compact 
manifold with and without boundary has been firstly derived in [11] and generalized in |12j . 
The high temperature limit for massless scalar fields in different topologies have been discussed 
by many authors [13 Mi 13 EH [23] since then. 

When concerning the Casimir contribution to the free energy (i.e. the Casimir free energy) 
in the case of material boundaries (e.g. parallel plates, pistons, rectangular box) all terms, which 
are of quantum origin like the Stefan-Boltzmann term are cancelled giving rise to the classical 
limit at high temperature [24" 1 l4"6 j [50] . 

Quite recently, the thermal Casimir effect was reconsidered for several types of fields in 
the static Einstein and closed Friedmann universe |51 |, I52 | . In order to obtain the Casimir free 
energy, it was proposed to use the renormalization scheme, which is usually applied in the 
case of material boundaries, also for the treatment of the topological Casimir effect at finite 
temperature. In fact, the Casimir free energy is the difference between the free energy of the 
topologically non-trivial manifold and the free energy of the tangential Euclidean/Minkowski 
space both filled with thermal radiation. Thereby not only the zero-temperature vacuum energy 
but even the finite-temperature contributions are renormalized. As a result the Casimir free 
energy tends to the classical limit at high temperatures, where the leading term is linear in 
temperature. 
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5 The equation of state 

Once the free energy T 8 is determined, the main thermodynamic functions can be computed 
by the following formulae: 

(5-1) 

Here II s is the internal energy, S 8 denotes the entropy and P e is the pressure of the photon 
gas. Since the thermodynamic functions are periodic under translations 9 4 6 + rh, for any 
rh £ Z bl( - X \ we can replace 9 by (6). 

Now we want to study the behavior of the free energy under a constant scale transformation 
/3 i — > A/3 and V i— > X n V with A £ R. Equivalently, this can be regarded as scale transformation 
of the metric g i— > X 2 g. A direct calculation yields 

vol x i ^ A%o/ x , 4V?) -> A- 2 ^(A?) 

^ A- 2 4\ ,f (Af ) A" 2 ^ (Af ) (5.2) 

with r, s = 0, 1. As a consequence, the free energy displays the following transformation behavior 



X 



imd^ )] 



V 





w 




(oi^/^V 1 ) 


G b2(x) (o^y^/^ 




bi(X) 










(5.3) 

showing that the free energy does no longer transform homogenously of degree —1. In fact, this 
shows that the free energy is not an extensive quantity. The violation of the scale invariance 
is caused on the one hand by the introduction of a scale ambiguity in the regularization of 
the vacuum energy and on the other hand by the harmonic modes of the Maxwell field. This 
anomaly has an implication for the equation of state. In fact, taking the derivative of (|5.3|) by 
A, setting A = 1 and using (|5.ip yields 



J*(P, V) = nP 9 ((3, V)V - /rV(/3, V) + T 9 (f3, V), 



(5.4) 



with 



r e (f3;V) = Res s= _ 1 [((s;A? 



n — 3 d 
+ ^5A |A=lln 



- £ (^)7Mm> 



j,k=i 



(5.5) 



6 fe2(x) (0|2Avr, 
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Together with the general relation T (/?, V) 
modified equation of state 



17° (0, V) - /3 _1 S (/?, V) one gets the following 



i 



U e (f3,V)-T e (f3,V) 



(5.6) 



Thus we have shown that the equation of state of a photon gas on a topologically non-trivial 
and compact spatial manifold differs from the conventional one in the large volume limit by the 
anomalous term T s . Irrespective of the topology of X, the "topological" terms in F e vanishes 
for n = 3. If in addition the zeta function is finite for s = — 5, then we obtain the familiar 
relation P e V = ^U 9 in three dimensions. 



6 Explicit results for the n-torus T n 

In the remainder of this paper we want to consider the photon gas confined to a n-dimensional 
torus X = T n in more detail. The n-torus is equipped with a flat metric 7 = Yli=x Lfdt 1 <g) dt l , 
where the sequence (t' l )f =1 denotes the local coordinates of T n and Li is the length in the i-th 
direction, so that the spatial volume is V = 11^=1 This configuration can be equivalently 
realized as field theory in a n-dimensional rectangular box X = [0, L\\ x . . . x [0, L n ] subjected 
to periodic spatial boundary conditions. 

Massless scalar fields in a box with periodic boundary conditions at finite temperature 
were studied in [131 C3 EH EH] based on Epstein zeta function regularization and in [22] using 
a multidimensional cut-off. An analysis of the finite-size effects in a universe with toroidal 
topology was presented in [IB] . 

We would like to stress once again that in the case of a non-vanishing boundary (where the 
fields in the rectangular box are satisfying Dirichlet or von Neumann boundary conditions) one 
has to take into account the contribution both from the inside of the cavity as well as from the 
outside region [24] for the calculation of the Casimir free energy. 

A Betti basis for U\{T n ) is given by pf ] = ±-i*vol^ n = dP, where if. T 1 T n is the 
canonical inclusion of the j-th position and j = 1, . . . ,n. Correspondingly, the dual Betti basis 

for H'%~ 1 (T n ) is provided by pf~ 1] = (-ly^p^ A . . . , Apf ] A ...,Ap£\ where the caret 

denotes omission of the respective element. By construction f Tn pp f\p£ ^ = 5jk- Furthermore 

pfk = L Lfc ijVoljn A i* k volj n = dP A dt k induces a Betti basis for %|(T n ), where 7, k = 1, . . . , n 
and j < k. As a consequence, W(Y n ) (p = 1,2) admit the following metrics 



hj} = diag 
h^l = diag 



V V 



15? = diaa I 



v v v - (6 ' 1 ' 



r2r2' r2r2 " ' ' r2 j2 



In order to obtain the free energy we have to calculate the spectrum and the zeta-function of 
A-i \imd*- It is easy to prove that the eigenvalues of A{ \i m d* are the sequence of real numbers 

^(A^" \i m d*) = Y^j=i {^~7~^ 1 parametrized by m G Zq := Z n \0. For each m G Zq, we 
introduce vielbeins e r j^, which are labelled by r and have components indexed by j, satisfying 
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n —In 

Ef 'jpn | . \ 2 ZTTTYlj . , . 

^rj,m ^r'j,m O rr ', ^nj,rn (^my^l limd^) J ~j~ j Ji ^* — 1, ... ,71. (6-2) 

j=l J 



Let us define the inner product ((p, %) 7 :=< </?, x >7 for 95, x £ Sl p (T n ,C), where x is the 
complex conjugate of Xi then the following 1-forms 



1 

n / t2\ 2 



: = E e ^rn ( ^ ] e 27r v /r TE?=i m iti ^ 7n G Zq , r = 1, . . . , n (6.3) 



w r n , . 



satisfy (ro r ^,ro r 'y) 7 = <5rT"5m,m' • A direct calculation gives 



, Jro r ^ if r= l,...,n-l 

d 2 G 2 diw r ^ = 4 (6.4) 
U it r = n, 



and 



Ai"|i mt i| = Vfh(Ai n \imd%) ™r,fh, mGZjJ, r=l,...,n-l, (6.5) 

proving that the set {w rr ^\r = 1, ... ,n— 1, mGZJ} provides an orthonormal basis of eigen- 
forms of Ai™|j m d* with respect to the inner product (,) 7 . Evidently, the multiplicity of each 
eigenvalue is just (n — 1). This is precisely the number of independent polarization states of the 
electromagnetic field. In terms of the Epstein zeta function (jB.ip we see that 

C(s;Af| w .) = (n-1) E n (s;^,...,^). (6.6) 
The regularized vacuum energy for the transverse modes (|4.22p is then given by 



£ reg * y 

'-vac • o „ 
2 T?— >0 



en\- 2 "n . 1 2tt 2ir /e/x\2»? 1 27r 2ir 

Y) En{ -2- 11] T 1 --r n ) + KY) ^ ( -2 +r?; ZI'-'V 



Applying the reflection formula (|B.5P and performing the limit rj — > leads to 



(6.7) 
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_n±l r , n + 1 n + 1 
2 r(— — )E n (— — ;Li, 



1 ^n) 



(6.8) 



By construction the vacuum free energy is finite which is - of course - confirmed by the fact 
that E n (^-; L\, . . . , L n ) converges and is positive for every n [53]. Moreover, for any given s, 
the minima of E n (s; Li, . . . , L n ) with fixed volume V = YYi=i appears at L\ = ■ ■ ■ = L n . 

As a consequence, there is no ambiguity in defining a finite vacuum energy. Inserting this 
result into (|4.26p . one immediately obtains the expression for the free energy on the n-torus, 
namely 
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n — 1 
n-1 



V vr-^r(^±i)£ n (^±i Ll ,...,L n) + ±± L){6, 



3=1 







E ln 



(fel,...,fcn) 



1 - e 



-2^ A /Er=i(^) 2 



!T ^i>i^4)-^ Mei ( o|W ^ /3 

/ j,k=l \ 

j<k 



1 



2tta 



V 



(6.£ 



Notice that TorH 2 (T n ;Z) = 0. For n = 1 we recover the free energy of pure electrodynamics 
on the circle (see e.g. [54]). 

For sake of completeness we want to sketch an alternative calculation of the free-energy 



Ti xT" 



which derives directly from (|3.35p . The eigenvalues of Aj/ |%p(ti xT") 1 are given by 



2 n /n , \ 2 



(o,o) _ ( 2irk \ sr^ f 2nk 

P ) ti\L 



r-r^m -Em , * 



an/ x 2 ( 6 - 10 ) 



1=1 



where = (A^> fci, . . . , k n ) £ Zq +1 := Z n+1 \0. Given the multiplicities of the eigenvalues, the 
zeta functions of the Laplace operators become proportional to the Epstein zeta functions in 
n + 1 dimensions, namely 



Ti xT" 



(n + lj^^g g) ifp=l 



By applying the duality relation ()A.3[) to 0i(^ |/3~ 1 ..) in (|3.35|) and using the Chowla-Selberg 
formula (|B.9|) by setting m = n + 1, I = 1 and ci = C2 = f^r--> c «+i = one immediately 



obtains (|6UD . 

In the torus case the vacuum energy is unique and the equation of state is modified only by 
the Riemann Theta functions. In contrast the free energy of the scalar massless gas transforms 
properly under scale transformations |23j. 

In the next two subsections we will derive explicit formulae for the thermodynamic functions 
in the low- and high temperature regimes, respectively. In order to simplify the calculations we 
will consider hereafter a uniform n-torus, i.e. L\ = ••• = L n = V~ . 



6.1 The low temperature regime 

It is evident that expression (|6,9h is suitable for studying the low temperature behavior (,8 > 1) 
of the system. Let us indicate this fact by a subscript and write J-f^iP; V) when the free energy 
is displayed in the form (|6.9p . 
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We can give an explicit series expansion for the last two terms in (j6.9|) . denoted by jf ow for 
further reference. Since | (#■,-) | < ^, the constraint (lA.5j) is satisfied. Eq. (|A.6j) can be applied to 
flow gi vlI1 g the following series expansion 



~P ^ ™ i_ e -^m (6-12) 

_ n(n- 1) y> 1_ l + 2(-l) m e L ^/ 3y " m 
2 ^ "i™ i _ e (2 w )2/3y^m 

where A is the number of components #j which satisfy \(9j)\ = \. Clearly < A < n. Apart 
from the first term, all others terms decrease exponentially for /3 — > oo. 

The vacuum energy ^Fn(V) is determined by taking the zero temperature limit of (|6.9p 



J*(V) :=\mxTf mu {p-V) = 

n — 1 n+i „,n+l s „,n + l , Tr i ,. 3-n (#i) 2 

= 2~ vr-— r(^-)£^(— ; 1, • • • , 1) V~ + V— (6.13) 

i=i 

If E Z n , the free energy of the Maxwell field equals the sum of the free energy of (n — 1) 
massless scalar fields (see e.g. the discussion of the massless scalar field with periodic boundary 
conditions in [17j and [23]). It is evident, that F^{X n V) = A -1 J#(V) for 6 e Z n . 
The internal energy derives from (j6.9|) and reads for the uniform n-torus 



a , , fl — n) 1 n+l ,71 + 1, ,71 +1 „ 2-n x— ^ (#1) 

tO/?; =^-^ F-«7r-— r(-±-)i^(^-; i,..,D+^Et 



+ E 2,(»-1)F-, W 9 ( 



(6.14) 



where we have defined the abbreviations k = (k\, . . . ,k n ) and \k\ = (X^ILi ^i) 1 ^ 2 ■ A direct calcu- 
lation using (|6.12p shows that the contributions in the last term of (|6.14p decrease exponentially 
for low temperatures. One obtains 

U°(V) := lim Uf^tf; V) = lim V) = J$(V). (6.15) 

/3->00 p— S>00 

The entropy of the system is given by 
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sL(P;V) =(i-n) J2 ln 



fcez" 



1 - e 



-2w/3V-n\k\ 



+ (n-l)V~njS ^ 



2vr|fc| 



fcez;? e (6.16) 



Using the series expansion (|6.12p and performing the zero temperature limit only the first term 
in (|6.12p survives, leading to 



lim S° ow ({3;V) = ln2 A . 

0— >OD 



3.17) 



This is a consequence of the fact that the ground-state of the system is degenerate of degree 
2 A . To verify this we notice that for a fixed the lowest energy level of n harm (see (|4.12|) ) in a 
given topological (monopole) sector reads 



n(n — 1) 
2 



'0,m 2 



(6.18) 



k=l 



n(n — 1) 

Here m = (mi, . . . , m^, . . .) G Z 2 labels the topological sectors. The degeneracy appears 
exactly for the parameter values 9j = 2 J » +1 , fcj € Z (i.e. = ^). Hence we have explicitly 

proved that the 3rd law of thermodynamics (Nernst theorem) holds for the photon gas on the 
n-torus. 

The pressure of the photon gas admits the following form: 



pL(P;V) 



n + l 

v — ~ 

n 



E 



2ir(n - 1)|A;| n — l _n±i n + l n + l 
5—1 tt~ n 2 1 v )^n{ „ ; 1 

2n/3V~n\k\ _ 1 2 2 2 



5 ? " " " 3 



(6.19) 



In the zero temperature limit, the pressure converges to 



P°{V) := lim PL(/5;^) = 

n — 1 , n+i n+i ^,n + 1. „ ,n + 1 , n — 2 aq-n) (6L-) 2 

For 1, one gets for the large volume limit limy_ s>00 Pq (V) = 0. Additionally, the following 
four cases can be easily distinguished: 

a) n = 1: Pq (V) = — is non vanishing whenever # ^ Z. 
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b) n = 2: The pressure is exclusively determined by the (vacuum) transverse modes of the 

gauge fields. Hence Pq(V) < for all V. 

c) n = 3: For 6 e Z 3 one gets P${V) < for all V. 

d) n > 3: For each 9 ^ Z™ there exists a critical volume T4 r it such that Pq (V cr it) = 0. For 

V < V cr it one has Pq{V) > whereas for V > V crit one has Pfj(V) < 0. 

A direct comparison of (16. 15ft and (16.20P shows that 

KiV) = ±Vg(V) + ^V^±'^f. (6.21) 
Alternatively this result may be obtained as the zero-temperature limit of the equation of state 

6.2 The high temperature regime 

An appropriate expression - hereafter called J r ^ i h (^;V) - for analyzing the high-temperature 
regime (/? <C 1) can be obtained from (I3.35P and (I6.1ip as follows: Taking m = n + 1, 2 = n, 
ci = c n = and c n+ i = ^ in ()B.9|) and applying ()A.3P one obtains after a lengthy 

calculation 



J* (/?; =(1 - n)^-^r(^±i)Ci ? (n + 1)/H" +1 V + ^ In/? + I InF + k 3 ln2vr) 



/3 



/3 



k. 



''» i-l = 1 fcGZg 



X if* 



2?rA;. 



n+l 



+ 



1 - 



2/3 



-K(0; 



Ln 



« E ln9 M°i : 



/3 



i,fc=i 
i<fc 



2vrF 



2tta 



5=1 



(6.22) 



n2 ~J n+8 and k 3 



n 2 -7n+4 



. These constants are 



with the three constants k± = - ~^ n+4 , «2 - j ... ( 
non vanishing for all dimensions n. 

We want to emphasize that the expressions for both Tf ow and J 7 ® igh are valid for all (3 and 
V, that is J~^ igh = J~f ow . However, their difference lies in the different polynomial structure as 
function of the temperature. 

Evidently, ()A.5p is satisfied so that for the uniform n-torus the sum of the last two terms in 
(|6.22p . denoted by /^ 5 / l (^; V), admits the following series expansion 
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+ 2 EE 



(27TT 

e 2 



m/3 1 V " 



cos(2irm(9j)) 



j=l m=l 



rn 



n — 2 

I — e (27r) 2 m/3- 1 V^ _ 



. (6.23) 



All terms show an exponential decrease for (3 — > 0. The same is true for the various derivatives of 
fhigh which appear in the different thermodynamic functions. Since K v (z) ~ \/j^ e ~ z (l + ^(i)) 
for \z\ — > oo [43J, the summands in the fourth term of (16. 22ft show an exponential decrease as 
P ->■ 0. 

As has been stated, (16.22P is exact. But what would be the result, if the general formula 
(|4.34p was used instead of (I6.22|) for determing the high-temperature limit? Taking into account 
(|6.1ip and the following explicit expressions for the Seeley coefficients of the Laplace operator 



Om(Af 



imd\ ) 



'(47r)~5 (n-l)V if m = 
1 — n if m = n 

if m/0,m^n, 



(6.24) 



we would have obtained (|6.22p up to the exponentially decreasing fourth term. 
In the high temperature regime the internal energy is given by 



U 9 high W; V) =n(n - l)vr-^ T<?±±) Cn(n + 1)/T (n+1) V + f 

n+4 n + 2 J^-v ^ — -\ H+1 , ->, 2-n / 2nV n , ->\ 

-4n(n-l)[3- — V^ £ £ \k\— Kn_ J -^-& n+1 1*| J (6 . 25) 

fc n+ l = l fc gZ n \ / 

where the relation 4-K v {z) = —K u ^\(z) — ^K u (z) have been used (see |43j). Apart from the 
leading Stefan-Boltzmann term all other terms in (|6.25p show a exponential decrease in the high 
temperature limit. In addition there exists a term linear in the temperature which is positive 
for all dimensions n. A lengthy calculation finally gives the following expression for the entropy 
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shghiP; v) = (™ 2 - r(^-) Cdn + i)/r n k + ^(1 - in/3)+ 

+ 2(n-l)/3-'Vi f; Yl k ! + i\^Kn (^k n+1 \k\\ 

-Mn-i) r ^v^ ± E^^^-if^Wil^l) (6 - 26) 

fcn+l = 1 fcGZJ ^ / 

72 — 1 

+ — 2~ £j,(0;l,...,l) -K3ln27r-K 4 lnV r 
where K4 = K2 + — -■ The pressure in the high temperature regime finally reads 



-fcfczii^*^ f E *»« 1*1^^1-1 (6.27) 

Q 

The term linear in temperature and volume is caused by the harmonic part of the gauge field 
on T n . It vanishes for n = 2 and is positive for n = 3,4. 



7 Summary 

To summarize, we want to highlight the main results of this paper: We shed some light onto the 
relationship between the topology of the compact, closed and connected spatial manifold X and 
the thermodynamic properties of quantum Maxwell theory at finite temperature. This has been 
achieved by providing a rigorous determination of the free energy for this system. It has been 
shown that there are inequivalent quantum theories giving rise to #-vacua, which are related to 
the first Betti number of X. The vacuum energy is the sum of the energy of the transverse modes 
and the energy of the harmonic modes of the Maxwell field. The finite temperature excitations 
are affected by the topology of X unless both the first and second cohomology of X vanish. 

In the Hamilton approach the quantization has been performed directly on the physical phase 
space using an appropriate parametrization of that space. The minimal subtraction scheme was 
applied to regularize the vacuum energy of the transverse modes. 

The functional integral quantization method on the other hand requires some advanced 
steps in order to overcome the Gribov problem: We constructed an integrable and globally 
defined measure on the space of all thermal gauge fields and introduced a consistent reduction 
procedure to eliminate the gauge degrees of freedom. By expanding the zeta function of the 
Laplace operators on Ti x X in terms of the zeta function associated to the corresponding 
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Laplace operators on the spatial manifold X we succeeded in proving the equivalence with the 
Hamiltonian approach. 

The thermodynamics of a photon gas confined to a n-torus is explicitly elaborated as a 
special example. Exact expressions for the thermodynamic functions for both the low- and high 
temperature regimes are derived and the validity of the Nernst theorem is proved. 

In light of the ongoing interest in the thermodynamic structure of quantum fields on spaces 
with non-trivial topologies, it will be a subject of future work to extend the present analysis 
to manifolds with boundaries and to apply these results for studying further the topological 
aspects of the Casimir effect in space-time models with compactified extra dimensions. 
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A Appendix: Riemann Theta functions 

In this section we want to summarize the main facts regarding the Riemann Theta function. 
Let us now introduce the r-dimensional Riemann Theta function with characteristics a, b £ W 
by 



a 



(u\B) = ^ exp {n\^l{m + a) f • B ■ (m + a) + 27r\/ z I(m + a) 1 " • (u + b)}, (A.l) 



for a symmetric complex r x r dimensional square matrix B whose imaginary part is positive 
definite, u G C r and the superscript f denotes the transpose. It has the modular property 



a + k 
b + l 



(u\B), k,l€Z r . 



(A.2) 



We write G r (u\B) := 9, 



(u\B) for the Riemann Theta function with characteristics (0,0). 



It can be easily shown that Q r (—u\B) = Q r (u\B) which leads to r 







(0\B) = Or 



By using the Poisson sum formula one can derive the following important "duality" formula 



(0\B). 



@ r (u\B) = det 



-iB) 2 a 



B- 1 ). 



(A.3) 



In order to provide a suitable asymptotic expansion of the Riemann Theta function we recall 
the existence of the infinite product representation (Jacobi triple product) which reads 



oo 



0i(u|£) = Y[{1 ~ q 2k ){l + g 2fc - 1 e 2 "^ 3Tu )(l + g^V 2 "^ 1 "), (A.4) 
k=i 
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where q = e w ^^ B . Let us now take B 
satisfy the constraint 



-It, with r£l and r > 0. If r and u = Ru+\/—lQu 



t - 2|3ftt| > 0, 



(A.5) 



then one can apply the series expansion of the logarithm and the geometric series formula to 
(|A.4|) , A straightforward calculation finally gives 



x X) 2km 00 00 



i) 



m+l 



fc=l m=l 
/ j ml g27rrr7i / / 



fc=l m=l 

oo 
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m=l 



J g27rrm 



(A.6) 



B Appendix: Epstein zeta functions 

In this appendix the main results on Epstein zeta functions [371 [38] (e.g. see also [231 [55] ) 
are summarized. Let us introduce the homogeneous Epstein zeta function E m in the variables 
(ci,...,c m ) eK m by 



Em{Si C\ , . . . , C r , 



E 

(fci,...,fe m )6Zg l 



i=l 



s G 



(B.l) 



where := Z m \{0}. The sum is convergent for 3?(s) > y. The analytic continuation of the 
Epstein zeta function is regular on the whole complex s-plane up to a unique pole at s = y. 
Under a constant scale transformation a h-> Aq, A£l, the Epstein function transforms as 



E m (s; Aci, . . . , \c m ) = A 2s £ m (s;ci, . . . ,c m ). (B.2) 
Epstein zeta functions in different dimensions are related by the Chowla-Selberg formula, 



E m (s;ci, . . . ,c m ) =£j(s;ci, . . . ,q) + — - — ^-E m _i(s - -;q + i, . . . ,c„ 

n«=i c * r ( s ) 2 

r(s) 



where 
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2vr s 



rii=i ^ 



E E 

(fci,...,fc i )e4 (fc ;+1 ,...,fc m ))ez™ 



J f k<\2 
Ci ' 
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2s-l 
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x K s _; | 2vr, 



n i EO 2 E (^*) 2 



(B.3) 



(B.4) 
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Here K v {z) denotes the modified Bessel function of the second kind [33]. The function s i— > 
T m i(s;ci, . . . , c m ) is analytic on C. Using the Chowla-Selberg formula one can prove the so- 
called reflection formula 

m 

7T ^ TTh TYl I 1 

Tr~ s T(s)E m (s;ci,. . . ,Cm) = -=n t (tt - s ) E m{— ~ S) — ), (B.5) 

1 li=l Cj 2 2 Cl c m 

from which the following results follows 

E m (0;c 1 ,...,c m ) = E 1 (0;c 1 ) =2^(0) = -1. (B.6) 

It follows from (|B.2p and (|B.6|) that the derivative of E m transforms under constant scale 
transformations according to 

E' m (0;\ Cl ,...,\c m ) = 2\n\ + E' m (0; Cl ,...,c m ), (B.7) 

where Cr( s ) = Sfc=i k~ s denotes the Riemann zeta function. By recursion we can express the 
Epstein zeta function in terms of the Riemann zeta function as follows: 

2 m ~ 1 -n^Tis - ±) 
E m (s; c 1 ,...,c m ) = 2cf% R (2s) + — - ]T K 2 > ( R (2s - i) 

1 W i=1 c i+1 U j= i cj 



2 4 



, „ m-1 -, oo [v* 

+ ^E^ E E (b. 



t II}=i c j ( felj ... ifei ) 6Z i fc< +1 =i [^+ic i+ i] s ^ 



x lf g _ i I 27r/c i+1 c 



i+l c i+lA 



Since lim s ^o(^j^y) = for a well-behaved function f(s) and using the reflection for- 
mula one obtains for the derivative of the Epstein zeta function E' m (0;ci, . . . ,c m ) := 
J^|s=0-Em(s;ci, • • • ,Cm) in s = 

7r""Tr(— ) ml 1 

£^(0;ci, . . . ,c m ) = ^'( ; c i 5 •••,£;) + — ™ — — -Em-z(ir; , • • • , — ) + r m i(0;ci, . . . ,c m ). 

H i=1 Cj 2 c/ + i c m 

(B.9) 
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